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We present the first study of the cusp effect (the movement of resonant poles due to the proximity 
of multiparticle thresholds) caused by the creation of diquark-antidiquark pairs. The cusp profile 
for such states is obtained from constituent counting rules. We compare the effectiveness of diquark 
cusps in moving resonant poles with that from a phenomenological form commonly used for meson- 
pair creation, and find that mesons tend to be more effective at lower energies {e.g., the KK 
threshold), while diquarks tend to be more effective at the charm [X(3872)] scales and above. 
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I. INTRODUCTION 

One of the most important discoveries of the past 
year in hadron physics was the experimental confirma¬ 
tion by LHCb [T] of a new type of hadron joining the 
qq mesons and qqq baryons - the tetraquark (qqqq) res¬ 
onance Z+(4430). In fact, this state is just one of 
an ever-growing menagerie of hadrons now believed to 
be tetraquarks, all of which have been observed since 
the 2003 discovery by Belle [2] of an unusually narrow 
charmonium-like state, the AT(3872). However, until the 
discovery of charged states (now called Zc), one could 
not be certain that the X(3872) was not just an unusual 
conventional charmonium state or a hybrid ccg; and un¬ 
til the observation of a phase S increase by ^ radians 
in the complex scattering amplitude in which Zc is pro¬ 
duced [1], one could not be certain that the states are 
true resonances as opposed to, say, kinematical reflec¬ 
tions of t-channel exchanges. There now seems to be little 
doubt that the (4430) is a tetraquark resonance with 
= 1+ quantum numbers and valence quark structure 
ccdu. In addition, all indications suggest that the charge- 
zero X(3872) is a = 1+ ccqq state, where qq is some 
linear combination of uu and dd. 

But waiting for the final confirmation of the tetraquark 
label did not preclude early speculation on the com¬ 
position and structure of such states. The initial in¬ 
terpretation - and still the most prevalent one - is 
that the tetraquarks are di-meson molecules, bound to¬ 
gether by color van der Waals-type forces (see [5] for 
a recent review). This interpretation is suggested by 
the proximity of the mass of several of the states to 
the corresponding two-meson thresholds. For example, 
TOjy( 3872 ) — rriDto — mjjo = —0.11 ± 0.21 MeV. How¬ 
ever, several other tetraquark candidates lie just above 
such thresholds - clearly muddling the simple bound- 
state interpretation - and several others lie nowhere near 
any two-meson thresholds. Furthermore, the substan- 
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tial prompt production cross section for the X(3872) at 
colliders seems to be incompatible with the state being 
solely composed of loosely-bound meson pairs HIS] - 
even taking into account a substantial modification due 
to final-state interactions [SIT]. Another well-known in¬ 
terpretation is hadro-charmonium [5], in which an ordi¬ 
nary charmonium state lies at the core of a light-quark 
cloud, although the cohesiveness of such states, and the 
extent to which they mix with conventional charmonium 
states, is unclear. 

In this work, we employ yet another noteworthy in¬ 
terpretation for_the tetraquark states, that of a diquark- 
antidiquark (S-S) bound-state pair. Originally proposed 
for charmonium-like tetraquarks in Ref. [5|, the diquark 
picture has the advantage of possessing a much richer 
color dynamics (the diquarks necessarily being color non¬ 
singlets), but it also has the potential drawback of pro¬ 
ducing many more tetraquark states than are currently 
observed. However, if one assumes that the dominant 
interactions are due to spin-spin couplings within each 
diquark, one can obtain a rather satisfactory accounting 
of the presently known tetraquark candidates m- 

One may wonder why the component quarks in a 6- 

5 bound state do not immediately rearrange themselves 
into color-singlet qq pairs, returning one to the molecular 
picture. Indeed, the strength of the color force between 
quarks (or antiquarks) in SU(3)-color representations Ri 
and i ?2 coupling to a representation R is proportional to 
the combination of quadratic Casimirs given by (72(7?) — 
C 2 {Ri) — (72(i?2)- The only attractive channels are the 
qq singlet (i?i = 3, R 2 = 3, R = 1) and the qq antitriplet 
(i?i = i ?2 = 3, i? = 3), with the former being twice as 
strong as the latter. 

Rather than treating the tetraquark as a metastable 6- 

6 bound state, we proposed in Ref. m a new paradigm, 
in which the tetraquarks are the quantized modes of a 
color flux-tube stretched between a rapidly separating S 
(color-3) and S (color-3) pair. This picture naturally ex¬ 
plains why many, but not all, of the tetraquarks lie near 
hadron thresholds (energies at which the color string is al¬ 
lowed to break); for example, the state A(4632) lies only 
slightly higher than the lightest charmed-baryon A+A” 
threshold at 4573 MeV, and it decays dominantly to this 




2 


baryon pair, as is explained by the fragmentation of the 
flux-tube to a light qq pair. The widths of the states be¬ 
low this threshold are not especially large, because they 
hadronize by wave function overlaps with the mesons 
formed from quarks in S and antiquarks in the S, which 
in turn can achieve a substantial separation (over 1 fm) if 
sufficient initial kinetic energy is imparted to the system, 
as in a B-meson decay. Evidence for this large separation 
is apparent in the decay of the .Zc(4430), which greatly 
favors coupling to tp{2S) tt rather than J/ipir [12 - even 
though the charmonium states have the same = 1 
quantum numbers and there is smaller phase space for 
the heavier (and spatially much larger) 

This paper presents the first dynamical study of the 
5-5 tetraquark picture, using the well-known constituent 
counting rules |13H23j . which determine the fall-off scal¬ 
ing of hard exclusive processes in the high-energy regime. 
In essence, the counting rules predict that the cross sec¬ 
tions, or equivalently, the form factors for processes at 
high values of Mandelstam s at fixed 9 cm, fall off as a 
power of s directly determined by the total number (in¬ 
coming plus outgoing) of fundamental constituents par¬ 
ticipating in the hard subprocess. Here we are inter¬ 
ested in comparing the effects of 5-6 state production 
with that due to meson-meson states on the masses of 
tetraquark resonances coupled to them, using dispersion 
relation techniques. In another paper to appear |24j . 
we will directly discuss the phenomenological uses of the 
cross-section scalings themselves. 

The method to be used is the well-known cusp effect, 
in which the presence of thresholds for the opening of 
on-shell states coupled to resonances creates a modifica¬ 
tion to the self-energy function that tends to drag the 
bare resonant pole mass toward the threshold. The basic 
idea appears to have been known since the early 1960’s, 
but was first presented in its current form in the mid 
1990’s [25], and first applied to the new heavy exotics 
in 2008 [55]. Furthermore, recent calculations m show 
that these states cannot just be cusps (although ones 
in the BB system might be [28l [29]); the presence of 
real resonances is required. In this paper we compare 
the effectiveness of this dragging effect due to the pres¬ 
ence of both 5-5 and two-meson thresholds, and find first, 
that the potential amount of shifting of resonant poles 
decreases for heavier-quark systems {KK vs. DD* vs. 
BB*), and second, that the 5-5 states become more ef¬ 
fective at pole-dragging than two-meson states for the 
DD* and BB* thresholds associated with the new heavy 
exotic resonances. Since 5 is not a color singlet, using a 
5-5 threshold in a QCD dispersion relation (where “on 
shell” usually means not only that the particles are not 
virtual, but also asymptotically free) must be interpreted 
with some care; for the purposes of this calculation, we 
assert that the substantial i5-(5 separation advocated in 
m creates states that are in a sense “almost” free, and 
therefore possess an on-shell threshold. 

This paper is organized as follows: In Sec. [IT] we briefly 
review the meaning and origin of the cusp effect, and 


establish the mathematical formalism used for its im¬ 
plementation, with some details relegated to the Ap¬ 
pendix. In Sec. [T^ we address in greater detail the issue 
of whether diquark pairs truly produce physically mean¬ 
ingful thresholds. Section [TV] presents a brief overview 
of the constituent counting rules used to establish the 
large-energy scaling of the 5-6 form factor. In Sec. 
we describe the algebraic results for both the mesonic 
and diquark forms, and present our numerical results in 
Sec. IVIl In Sec. IVIII we summarize and indicate future 
directions. 


II. THE CUSP EFFECT 


The cusp effect has a rather straightforward origin in 
the analyticity of the self-energy functions n(s) that ap¬ 
pear in the propagator of resonant or bound states and 
source the creation of their decay products. Closely fol¬ 
lowing notation introduced in Refs. [25l[2^, we start with 
the propagator denominator 

P-^{s) = iMl^-s)6c.g-nc,g, ( 1 ) 

where a, /3 index resonances that can mix; for simplic¬ 
ity, in this work we assume only unmixed propagating 
states and henceforth suppress this index. The sign of 
n(s) is chosen to match that appearing in the majority 
of quantum field theory texts {e.g., [3D]); when positive, 
it is seen to correspond to an attractive interaction. 

n(s) is analytic everywhere in the complex s plane ex¬ 
cept for cuts (or poles) along the positive real s axis that 
result from the opening of on-shell channels, which we 
label by i. Consider explicitly the creation of such two- 
particle states of masses mi^i, m 2 ,i via the form factors 
Fi{s), conventionally normalized to unity at threshold, 
Sth.i = (wi_i -|-TO 2 .i)^- Defining the (dimensionful) cou¬ 
pling constant to state i as gi, one has 

Imn(s) = ^g'^i p^{s) Ff{s) e{s - sth,i) , (2) 


and Pi is the two-body phase-space factor given in terms 
of the c.m. momentum or the Kallen function A by 


~ ffs s 


(3) 


Defining Sth.i as the lowest threshold (and therefore the 
branch point whose cut extending to s = oo overlaps all 
others), one may apply Cauchy’s theorem to a contour 
that goes around this cut, obtaining the standard unsub¬ 
tracted dispersion relation: 


1 pOO T 

Ren(s) = -P/ ds' — 

^ Js,. 1 S 


, Imn(s') 
— s 


(4) 


where P indicates the standard Cauchy principal value 
prescription. For reasons to be discussed below, we find 
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it useful to use the once-subtracted form (at s = 0) of 
Eq. (|4|: 


Ren(s) 


Ren(O) + -P 

•'Sth, 1 


ds' Imn(s') 


(5) 


One notes immediately from Eqs. ([^-(|^ that Im n(s) 
is zero on the real s axis for s < Sth,i, and positive for 
s > sth,i- One further sees from Eq. Q or Eq. ^ that 
the integrand in Ren(s) is positive for all s < Sth,i since 
s' — s > 0 in the entire integration range, and there¬ 
fore acts as an attractive interaction. As shown below, 
explicit functional forms for F(s) give rise to a positive 
Ren(s) that reaches a peak, or cusp, at Sth,i and falls off 
in either direction, remaining always positive for s < Sth.i 
but generally passing through the axis (and effectively 
creating a potential barrier) for some value of s > Sth,i- 
This effective attraction acts to pull the pole position of 
the propagator Eq. Q toward Sth,i, a synchroniza tion of 
the resonance with a threshold (although in Sec. |VI| we 
see some interesting exceptions to this expectation). 

In the bulk of this work, for simplicity we specialize to 
the equal-mass case mi^i = m 2 ,i = rrii, which is rather 
closely satisfied for the cases of experimental interest. 
Full expressions analogous to the ones appearing in the 
text below are presented in the Appendix, where it is 
seen that the relevant expansion parameter is 


/ - 1712, 
\mi^i + m2,i) 


For A(3872), with mi^i = and TO 2 y = mjjto, one 
finds e = 1.35 • 10“^. In the equal-mass case. 


P^{s) = 




(7) 


and its analytic continuation for s < Sth.i = 4m| is 


^'^(s) = 




( 8 ) 


To illustrate the cusp explicitly, let us consider the 
special case Fi{s) = 1 described in [2S]. This particular 
choice is unphysical because unitarity and perturbative 
QCD require real hadronic form factors to fall off at large 
s, but it is an interesting toy model for cases in which the 
dispersion relation is dominated by the form factor near 
threshold. It is also interesting because the dispersion 
integrals can be carried out explicitly; the nonvanishing 
large-s behavior of Fi{s) requires one to use the once- 
subtracted form Eq. ([^ to obtain finite results at finite 
s. Allowing the constant Re 11^(0) = — :| to renormalize 
the pole mass Mq in the propagator Eq. Q, one finds: 


Reni(s) Pi I-Pi 2pi 1 

- 2 - = ~ In —-=-tanh Pi, s> Sth.i, 

gf TT 1+p, TT 


2% 


= -Vi + — tan Vi, s < Sth.i ■ 

TT 


(9) 



FIG. 1: The threshold cusp in Re Iii{s)/g‘i (green, solid) 
and corresponding (red, dashed) for form fac¬ 

tor Fi{s) = 1, as given by Eq. Q. s is expressed in units of 
s/sth,i, and /gi is shifted to equal 1 at s = Sth,i. 


Using Eqs. Q-Q, note that the tan“^ and tanh“^ 
terms are analytic continuations of the same function of 
s in different regimes; the only unmatched term is —Vi, 
which is induced by analyticity of n(s) and the discon¬ 
tinuity of Imni(sth,i)- In Fig. we plot Eq. ® [and 
Imni(sth,i)/5|] as a function of s/sth,i, and see that the 
below-threshold term —Vi is responsible for the cusp. 

Cusp behavior also appears for more realistic choices of 
form factor Fi{s), although the dispersion relation inte¬ 
grations must often be carried out numerically. Two fea¬ 
tures of the integrals hamper the convergence of numeri¬ 
cal calculations: First, the integration range stretches to 
s = oo, and second, the denominator factor s' — s pro¬ 
duces a logarithmic singularity when s > Sth,i- A simple 
change of variable cures the first problem; noting from 
Eq. 0 that s € oo) maps to pi € [0,1), one defines 
the integration variable p' analogously to Eq. 0, but 
with s —>■ s'. Then one finds the unsubtracted relation 
Eq. 0 to become 


Ren^(s) _ 2 sth,» 

gf TT s 




P? 


1 - P? P? 




( 10 ) 

where [inverting Eq. 0] s' = Sth,i/(1 — p'i^)- In fact, the 
once-subtracted relation Eq. 0 becomes even simpler: 


4 [Ren,(s) - Ren,(0)] = - p / dp) ^ Ff (s') . 

9i JO Pi Pi 

( 11 ) 

Since, as discussed above, one expects Reni(s) —)■ 0 as 
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s —> oo for physical -Fi(s), no information is lost by us¬ 
ing the once-subtracted and functionally simpler form 
Eq. (11). The logarithmic singularity does not arise for 
s < sth.i) since pf = —vf according to Eqs. (^-(|^, one 
may rewrite the integrand factor in Eq. (11) for s < Sth,i 
as the nonsingular form 


P? 


^'2 


^ 1 - 


V2 


,,2 ■ 


( 12 ) 


The logarithmic singularity for s > Sth.i is readily han¬ 
dled via an integration by parts on the l/(s' —s) factor in 
Eq.lU), followed by conversion to the integration variable 
p'. One obtains: 

^ [Reni(s) — Re 11^(0)] = 


dp'. 


Fi{s') + 


dF^s') 


{1-p'^y ds' 


In 


\P? 


Pi\ 


1-P. 


(13) 


The singularity at pi = remains, but now it is inte- 
grable for any smooth Ff{s) that asymptotes to zero as 
s —>■ oo. Moreover, the non-logarithmic part of this in¬ 
tegral is just the exact differential d[pi{s')Ff{s')\, whose 
argument vanishes at s' = Sth,i and also at s' = oo if 
F^{s') falls off faster than l/pi(s'); in such a case, the 
(p'-independent) term 1 — in the logarithm may also 
be deleted. Since the function obtained from Eq. ([^, 
i.e., the explicit left-hand side of 0 or ( |I^ , is cen¬ 
tral to our analysis, let us henceforth abbreviate it as 
TTiis) = [Ren,(s) - Reni(0)]/5f. 


III. DO DIQUARK PAIRS HAVE 
THRESHOLDS? 

Before proceeding, one must address the issue men¬ 
tioned in the Introduction of whether a S state, being 
colored, really possesses a well-defined mass, and there¬ 
fore whether the 6-S state really possesses a well-defined 
threshold. Strictly speaking, the quoted masses of col¬ 
ored particles like quarks depend upon the choice of 
renormalization scheme, and therefore do not carry the 
same status of being observables as do meson masses. 
Nevertheless, quark masses are extracted using a vari¬ 
ety of methods (quark models, lattice gauge theory, sum 
rules, etc.), and the numerical values thus obtained are 
bestowed with some degree of physical significance. 

An elegant discussion of such issues appears in the 
classic text by Georgi [31] (Sec. 3.2); he considers a uni¬ 
verse in which the strong coupling constant as is as small 
as oemi but confinement still occurs. The scale Aqcd 
then becomes extremely large (~ 10^° cm), and for all 
practical purposes, quarks would be just as visible in 
experiments as electrons, and would possess easily de¬ 
termined quantum numbers such as mass. Georgi then 
continues the analysis by increasing Aqqd towards its 


physical 0(250 MeV) value, and argues that not until 
Aqcd — asTUq (including the running value of a^) do 
quarks become inseparable from their hadrons, and hence 
hidden from view. 

We take this lesson as our touchstone, that substan¬ 
tial physical separation implies identifiability as particles 
with well-characterized masses, even in the presence of 
confinement. The color dynamics allows diquarks to form 
as bound states, and kinematics permits them to separate 
a distance before being forced to hadronize. The salient 
question then becomes: How far is far enough, before 
the diquarks can reasonably be said to appear as iden¬ 
tifiable particles possessing a well-defined mass thresh¬ 
old? We argue that this event occurs when the 6-S state 
can no longer be mistaken for one in which the diquark 
wave functions have substantial wave function overlap. 
Of course, the size of diquark wave functions is not at 
all a known quantity, but we can obtain a reasonable es¬ 
timate by considering the size of mesons with analogous 
quark content. That is, instead of cu and cd diquarks, we 
consider D mesons. Again, meson radii are not directly 
observable quantities and the precise nature of the fall-off 
of their wave functions with r is unknown, but at least in 
this case some calculations have been performed. For ex¬ 
ample, Ref. [32] calculates the electromagnetic charge ra¬ 
dius of the to be 0.43 fm, rather smaller than charge 
radii of the unflavored mesons, due to the presence of 
the heavy c quark. One may expect the diquarks to be 
slightly larger because the initial color attraction between 
the quark components is smaller; however, the true size 
is dominated by some combination of nonperturbative 
gluodynamics and the heavy quark mass, so we expect 
the diquarks to be not much larger. We take from this 
result that 1 fm is not an unreasonable estimate for the 
onset of significant separation between the diquarks, and 
hence the identification of a well-defined diquark mass 
and pair-production threshold. Again, an important in¬ 
gredient in this estimate appears to be the presence of 
heavy quarks, which might very well be related to the 
reason why exotics have first become clearly visible in 
the charm system. 

Furthermore, as first pointed out in Ref. [l^, the cusp 
effect treatment of the previous section holds at full 
strength [as in Eq. ([^] only for thresholds of particles 
in a relative s wave. Only certain 5-5 states will fall into 
this category, in precisely the same way that the con¬ 
ventional cusp effect is most prominent for thresholds of 
mesons produced in a relative s wave. 

A very interesting question is how the diquark thresh¬ 
olds disappear if their separation is insufficient to pre¬ 
vent a large wave function overlap. We speculate that 
these thresholds gradually dissolve into correlations be¬ 
tween the couplings to meson pairs to which the exotics 
preferentially decay, and such correlations would create 
a smeared “cusp” that would gradually disappear for 
smaller and smaller diquark separations. The exact sepa¬ 
ration at which such a structure first appears, and hence 
the assumption that a viable diquark cusp effect is sup- 
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portable, is very much open to debate, but we believe 
1 fm is a reasonable estimate. The diquark cusps stud¬ 
ied in this paper are of course based upon an idealized 
picture of this behavior. 


IV. CONSTITUENT COUNTING RULES 


In this section we present a brief summary of the con¬ 
stituent counting rules developed over a number of years 
following the advent of perturbative QCD (pQCD) [131- 
E5] . They reflect the underlying conformal features 
and scale invariance of the QCD coupling. More re¬ 
cently, they were derived nonperturbatively by Polchin- 
ski and Strassler [33] using AdS/QCD. A more thorough 
pedagogical introduction, particularly with an eye to¬ 
ward discerning exotic hadronic structure, is presented in 
Ref. [34]; here, we present an abridged discussion identi¬ 
fying the central points leading to the correct counting. 
In Ref. [3S] , the same authors show how to use the large- 
s scaling behavior to study the underlying generalized 
parton distributions and distribution amplitudes enter¬ 
ing into these processes. 

The constituent counting rules for large-angle scatter¬ 
ing processes (i.e., those in which no kinematical variable 
is small due to the near-collinearity of any particles) at 
high c.m. energy originate from a remarkably simple 
source: Each of the individual constituents being scat¬ 
tered must be redirected through a finite angle by a large 
momentum transfer. In this limit, all constituent masses 
are negligible, and all three Mandelstam variables s, t, 
and u ~ —s — t are large, so that one may express all 
dimensionful quantities in powers of s and dimensionless 
coefficients as functions of t/s. In the pQCD picture, 
these deflections are accomplished through hard gluon 
exchanges; if leptons are included, then hard electroweak 
boson exchanges also appear. In the case of AdS/QCD, 
the counting rules reflect the twist dimension of the in¬ 
terpolating field at short distances. 

For the moment, let us consider processes in which 
the constituents are all fermions. Then, assuming that 
each of the n = n-m + Uout constituents shares a finite 
fraction of the total s, the leading-order Feynman dia¬ 
grams require a minimum of ^ — 1 hard gauge boson 
exchanges and a minimum of ^ — 2 internal constituent 
propagators. These features supply factors of 

and 1/y/s^^ to the invariant amplitude Ai, respec¬ 
tively. Noting that each external constituent fermion 
field carries a spinor normalization scaling as s *, one 
sees that all of the fermion scaling factors cancel except 
for an overall factor s. In total, 

M (X . (14) 

Assuming a conventional scattering process in which the 
constituents combine into two initial and two final parti¬ 


cles, the cross section is given by 


da 

dt 


1 




1 

gn-2 


f 



(15) 


where / has the appropriate mass dimension (M^*^”®) 
to match that {M~^) of the left-hand side, but does not 
itself scale as a power of s; its dimensionful factors are 
essentially the amplitudes describing the binding of the 
constituents into the composite states, i.e., decay con¬ 
stants. 

Each external gauge boson introduced (e.g., turning an 
electroproduction process into photoproduction) removes 
two external fermion lines (-s/s)^] and one hard gauge 
boson propagator (~ 1/s), leaving the form of the scaling 
formula for A4 and da/dt invariant. 

Hadronic form factors in the large-s regime can be 
studied analogously. For example, according to Eq. (141, 


the electromagnetic (or any other current-produced) 
form factor Fx{s) appearing in the pair-production am¬ 
plitude M for tetraquark states X scales as 


Fx{s) 


1 


g^(l+l+i + i )-2 53 


(16) 


While the origin of 1/s factors via hard exchanges is nat¬ 
ural and simple, one may worry about technical compli¬ 
cations in real perturbative QCD that disrupt the sim¬ 
ple counting. These effects include the running of q:s(s), 
the renormalization-group scaling of the constituent dis¬ 
tribution amplitudes [361 El], the presence of Sudakov 
logarithms [55MD] . “pinch” singularities due to the van¬ 
ishing of internal gluon propagators im, and endpoint 
singularities occurring in configurations where some of 
the constituents do not share an 0(1) fraction of the hard 
scale s [H]. However, it is believed that the net result 
of these effects is not sufficiently severe as to change the 
leading s power scaling for exclusive processes. 


V. MESONIC VS. DIQUARK FORM FACTORS 

The discussion of the previous section shows that the 
pair-production form factor Fi{s) for a state with n con¬ 
stituents in the high-s regime scales as 1/. Of course, 
this scaling is not expected to hold in the low-s region, 
particularly since confinement physics is not taken into 
account in this fundamentally perturbative approach. 


A. Meson Form Factors 

In the case of pair creation of conventional mesons, 
an exponential form is traditionally favored by phe¬ 
nomenological fits to data. For example, for KK pro¬ 
duction, Ref. [26] uses Ff = exp(—fc/i?^/3), where 
R = 0.6 fm, while for BB* production. Ref. [29| uses 
Ff = exp(—s//3/) (once the coupling constant gi is 
removed), and studies the choices fdi = 0.4, 0.5, and 
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0.7 GeV. To put these choices on a common footing, note 
from Eq s. (|^ and Q that ki = — Sth,i, and that the 

form in |29) can be normalized to unity at threshold by 
multiplying by a constant, = exp[—(s — sth,i)//3z^]- 
Then the exponential form may be written as 


if (s) = exp 



where 


tJ-i = 


12 


w 


(17) 


(18) 


so that fii = 0.738 for [26], and R — 0.6 fm corresponds to 
Pi = 1.14 GeV. If ki from Eq. (|^ is identihed with the 
nonrelativistic momentum |ki|, then one may compute 
the normalized nonrelativistic matter density p(r) (not 
to be confused with pi) as its Fourier transform: 


giving in this case a Gaussian form: 


(19) 


( 20 ) 


Inserting Eq. 0 into Eq. ( |TI| ) leads to integrals for 
Trfs) that do not appear to be expressible in closed form 
except at special values of s. That the subtraction is 
performed at s = 0, and hence that TrfO) = 0, is guar¬ 
anteed by noting that s = 0 in Eq. ([^ corresponds to 
vf = —Pi = oo. In terms of the Kummer (confluent hy¬ 
pergeometric) function U{a,b,z), the function Trfs) and 
its slope 7r'(s), can be computed at a few key values of s: 


7rf0) = 0, 

7J’i(sth,i) = 0, p) , 

Trfoo) = -^C/(i,l,/r), 


7r'(oo) = 0. 

( 21 ) 


These values can also be written as somewhat more com¬ 
plicated expressions in terms of the modified Bessel func¬ 
tions Kq and Ki. Note especially the cusp 7r'(s((, ■) = oo 
(approaching from s < Sth,i), which follows directly from 
Eq. ( [ll| ) for any F^{s) that is smooth at s = Sth,i (Pi = 
0). Tlie slope 7r'(sf J (approaching from s > Sth.i), on 
the other hand, is generally finite (and negative) due to 
the smoothing effect of the principal-value prescription 
in Eq. ( |TT| ). 

It is important to note that a form factor F^{s) ex¬ 
ponential in s cannot truly represent the full physical 
amplitude |43j in the entire complex s plane, since it pro¬ 
duces an essential singularity for large s in some direc¬ 
tions. Such behavior in dispersion relations would lead, 
for example, to a violation of causality. For our purposes, 
however, one may suppose that the behavior of Ef (s) re¬ 
mains numerically close to the exponential form along a 


substantial portion of the real s axis, but that the exact 
form contains a functional dependence in s (for example, 
the power-law fall-off predicted by quark-counting rules) 
that restores the proper analytic behavior for all complex 
values of s. 


B. Diquark Form Factors 


The discussion of Sec. m shows that the tetraquark 
form factor F{s) at large s scales as 1/s^ due to the ex¬ 
change of hard gluons needed to maintain the integrity 
of the exclusive tetraquark state. Indeed, the true scal¬ 
ing is [as(s)/s]^; however, the large-s scaling of Us is 
logarithmic and therefore varies slowly compared to the 
power-law behavior, and henceforth is neglected. 

The corresponding large-s behavior was not used for 
the mesonic case because one expects the four quarks pro¬ 
duced at any low or intermediate s immediately to con¬ 
fine into two hadrons; the only color interactions between 
the two mesons are then the “color van der Waals” forces 
represented by final-state interactions. In the model of 
Ref. m, however, the (colored) diquarks separate a sig¬ 
nificant distance before they are forced to hadronize, and 
so the fundamental color forces can be expected to re¬ 
main active at much lower values of s than for the meson 
case. Since the natural scale at which the high-s scaling 
should become significant is given by the diquark mass 
^Sth,i = 2 to 5 , we model the diquark form factor by 

E,(s)=(^)', (22) 


which of course is properly normalized, Ei(sth,i) = 1- 
Using Eqs. (§ and ( [T^ , the matter density associated 
with this form factor is 


where rc is the diquark Compton wavelength llms- 
In this case, the integrals in 'iTi{s) can be performed in 
closed form. First, for s > Sth.i, one hnds 





1 + Pi A 

1-pJ 

^ 4_ 
or Pi 


2048 

3003 



(24) 


In the s < Sth.i case, we define [using Eq. ([^] the param¬ 
eter 


s 

li = - 

5th,i 


1 


(25) 


which varies from 0 —> I as s increases from 0 —>■ Sth,i- 

















7 


One then computes: 


7ri(s) = — 6 
7r7° 


1 - 7z 
7i 


,-i 


5-7* 


— — tan , ,_ 

2 V V7i(l - 7,) 


o , 1 2 2 8 3 16 4 128 5 

^ I 3^* 15^* 105^* 315^* 3465^* 


256 

9009 


7 i 


(26) 


These precise forms are not terribly illuminating in their 
own right, but they can be used to compute various limits 
for purposes of comparison with the mesonic case. For 
example, Eq. (26) appears to have a very strong singu¬ 
larity at 7 i = (l(s = 0 ), but the complicated polynomial 
in 7 i has precisely the right form to cancel the all terms 
in the expansion of the tan“^ expression in powers of 7 ^ 
to the same order to give 7ri(0) = 0, and furthermore to 
give a finite value at 7 i = 1 (s = Sth,*) that matches the 
one from Eq. (24) at pi = 0, but with a derivative from 


below that is infinite. In comparison with Eq. (21), 


7r*(0) = 0, 

7’‘i(sth.i) = 300371-■ 
^i(oo) = —900977 


2048 


Tr'fOl = __ 

45 0 4571-sth,, 

<(St"h,7) = 00 > 

7r'(oo) = 0. 


(27) 


VI. NUMERICAL RESULTS 

In this section we present explicit numerical compar¬ 
isons of the effects of threshold cusps due to mesons com¬ 
pared to those due to S-S pairs. But first, we make a few 
comments about how best to present and interpret these 
effects. 

While it is rather suggestive to describe the attractive 
cusp in Reni(s) at s = Sth,i as a potential well that at¬ 
tracts a bare pole to the threshold, it is important not to 
be overly seduced by the analogy with a configuration- 
space potential well, which always attracts a particle to 
the vicinity of its minimum. The intuition one develops 
for the cusp effect needs to be more nuanced. In fact, 
the entire portion of the cusp over which Re 11 ^( 5 ) is pos¬ 
itive and appreciable in magnitude provides a source of 
attraction for the pole. For example, even if the bare 
pole mass sits exactly at Mg = y'Sth.i, the cusp can drive 
the pole mass to a value slightly larger than such 

an effect is visible in the results in Table II of Ref. [2^ . 

The corrected pole mass Mpoie is the pole of the propa¬ 
gator, be., the zero of Eq. 0- As the physical resonance 
pole first appears on the second Riemann sheet, one must 
take care to choose the sign corresponding to the correct 
branch in the self-energy function 11 ^( 5 ). The effect of 
the threshold cusp on the position of a pole can then be 
described succinctly in terms of the dimensionless ratios 
Mo/ysth,i and Mpoie/\/sth,i- In light of the observation 
that the whole profile of the cusp is active in dragging 
Mq to Mpoie, one expects that the cusp is most effec¬ 
tive in dragging the pole if either the coupling gi to the 


opening threshold is large, or the cusp function Re 11^ (s) 
is broad in s. The question of whether the mesonic or 
diquark threshold cusp is more effective in pole dragging 
thus comes down to these two criteria. 

To exhibit the effectiveness of the cusp numerically, 
we first superimpose two plots in Fig. The first is a 
sample (diquark) cusp 7ri(s) with gi = 0.370 GeV and 
■<ysth,i = 3.872 GeV, which uses the abscissa x = s/sth.i- 
Its peak is normalized at threshold {x = 1) to match that 
of the second plot, which uses the abscissa x = 
and an ordinate that gives a measure of the effectiveness 
of the cusp in dragging the pole, chosen to be 


Mpoie - Mq 
Sth,* 


(28) 


The pole-dragging function y very closely follows the 
shape of the cusp, but does not precisely match it. Qual¬ 
itatively, this result means that the cusp does indeed at¬ 
tract the pole over a significant range of Mo/ysth.*- To 
consider the effect more carefully, we characterize regions 
of the plot: 

^ Mq ^ Mpoie ^ 

^ -Sth,* 

attraction of below-threshold pole toward Sth,*; 

2 Mq Mpoie . 

^ Sth,* ^ Sth,* 

attraction of below-threshold pole past Sth,*; 

3 1 < Mq ^ Mpoie . 

^ Sth,* ^ Sth,2 

repulsion of above-threshold pole from Sth.b 

^ ^ ^ Mpoie ^ Mq ^ 

Sth,* SthjZ 

attraction of above-threshold pole toward Sth,*- 

The other two possibilities, be., repulsion of a below- 
threshold pole from Sjh ^ and attraction of an above¬ 
threshold pole past Sth.i, do not occur for cusp functions 
of the types considered here. The interesting possibility 
(Region 2) mentioned above, of the cusp causing a below- 
threshold pole to overshoot threshold, is represented in 
Fig.i as the sliver of the plot between the line \ — x and 
the vertical line x = Mo/.ysth,i = 1. At these input val¬ 
ues, the peak of the pole-dragging function y lies slightly 
below the peak of the cusp function (at a: = 1), but still 
above the dividing line 1— a;, so that Region 2 has a fi¬ 
nite extent. Equally interesting is Region 3, between the 
lines a: = 1 and where y = 0, in which the bare pole 
Mq lies near threshold but is repelled by it nonetheless. 
Both of these effects are caused by the continuing attrac¬ 
tion (to larger s values) of the cusp beyond Sth*- From 
Fig.[2j one sees that the most effective dragging of the 
pole (the maximum of y) falls very close to the dividing 
line y = 1 — a;, which is the point where Mpoie = y^sth7- 











































0.0010 


0.0010 



FIG. 2: A sample threshold cusp (solid, green) using the di¬ 
quark form presented in Sec. |VB1 plotted as a function of 
X = ^Js/ Its peak normalization is chosen to match 
that of the overlaid (dashed, red) pole-dragging effect of 
this cusp, y = (Mpoie - Mo)/^Sth,i, with gi = 0.370 GeV 
and y/Sth,i = 3.872 GeV, and plotted as a function of a 
second abscissa x = Mo Also plotted is the line 

y = l — x=l — Mo /y^Sth,i (dotted, blue) dividing Regions 1 
and 2, as defined in the text. 


FIG. 3: Comparison of the effectiveness of pole dragging by 
cusps [y of Eq. (281] as a function of a; = Mo /^Sth,i, as created 
by diquark (solid, green) and meson (dashed, red) thresholds. 
The values of (in GeV) corresponding to the mesonic plots, 
in increasing width of the prohle (or darkness of the shading), 
are 0.5, 0.8, 1.1, 1.4, and 1.7. The mesonic plots are scaled to 
the same cusp height, in the manner discussed in the text. 


We now turn to a direct comparison between the pole¬ 
dragging efficiency of diquark and mesonic cusps, which 
are chosen to have the same value of Sth.i- Of course, 
the presence of a cusplike structure not coinciding with a 
known mesonic threshold could suggest the presence of a 
distinct diquark threshold, possibly mixed with a S-S res¬ 
onance as discussed here, or even a heretofore unknown 
qq resonance; and even in the case that a meson and di¬ 
quark threshold coincide or overlap, the shape of their 
combined effect would be distinct from that of a meson 
threshold alone. For the diquark form, the only free pa¬ 
rameters are the coupling constant gi and the bare pole 
mass as a multiple of threshold, Mq /y/Sth,i- The mesonic 
cusp, on the other hand, contains the additional free pa¬ 
rameter Pi indicative of a typical hadronic scale, which 
should assume roughly the same « 0.5-1.0 GeV value 
for any hadronic system; alternately, it can be expressed 
[Eq. (18l] as the dimensionless combination gi, which 
changes from system to system depending on the value 
of Sth.i- 

In Fig. we compare the pole-dragging effectiveness 
parameter y of Eq. (28) as a function of x = Mq /^Sth.i 
for the choice ^ysth,i = muo+mjyto = 3.872 GeV relevant 
to the X(3872) . The diquark pole-dragging plot is identi¬ 
cal to that in Fig.(it is the dashed curve there), where 


it corresponds to diquark = 0.370 GeV. The mesonic 
pole-dragging plots are presented for several values of Pi. 
Since, as discussed below, the values of gi,diquark are cho¬ 
sen to scale with gi,meson so that the diquark and mesonic 
cusp functions have the same height in each case, and 
since the absolute height of each mesonic pole-dragging 
curve changes with Pi, one should really view Fig. as a 
family of plots: a diquark and a mesonic curve for each 
value of Pi- The diquark plots inherit Pi dependence only 
through this scaling. The diquark and mesonic plots give 
the indicated value of y only for /3i = 1 GeV (where 
gi,meson = 0.474 GeV), and both would be proportion- 
aliy larger (smaller) for Pi > 1 GeV (< 1 GeV). How¬ 
ever, for fixed ^Sth,i the diquark curves have a fixed 
shape, and hence are scaled to the single (solid) curve 
in Fig. 1^ while the mesonic curves become broader as Pi 
grows. They increase monotonically in width (and hence 
in pole-dragging effectiveness) as a function of Pi, achiev¬ 
ing near-parity with the diquark curve at Pi = 1.7 GeV. 
However, such a large value seems inconsistent with that 
expected from ordinary hadronic matter; if one limits to, 
e.g., values Pi = 1.1 GeV or less, then the mesonic curves 
are much narrower than the diquark ones scaled to the 
same height, and hence mesonic thresholds do not drag 
resonant poles as effectively as diquark thresholds. 

However, this very interesting result is also dependent 
upon the absolute size Sth,i of threshold. One expects 
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the same range of /3i values to occur due to hadronic 
confinement physics, but the mesonic plots depend upon 
the dimensionless combination of Eq. ( |18| ). In the 
case of the KK threshold studied in Ref. |26j ! ^/sth^ = 
0.991 GeV, while for the studies of Zi, states in Ref. [53], 
the threshold is .^Sth.i = itib + ms- = 10.604 GeV. 

One must also take into account the chosen size for the 
coupling constant gi. In Fig. [^ the same cusp function 
peak value is used for both diquark and meson forms. Of 
course, no diquark pair production coupling constant has 
ever been measured experimentally; however, as noted in 
the Introduction, the fundamental strength of the color 
interaction between a qq pair forming a 3 is fully half 
as strong as that between a qq pair forming a 1, so it is 
reasonable to take the diquark and meson gi values to 
be comparable in size. In the plot captions, gi,meson is 
the given numerical value, and gi,diquark is derived from 
normalizing the peak of its cusp function to match that of 
the mesonic form; this ratio is computed using the values 
of 7ri(sth,i) in Eqs. (211 and (27). 

Furthermore, the coupling gi corresponds to the cre¬ 
ation of hadrons containing heavy quarks of species Q 
{s for KK, c for DD*, b for BB*), and therefore is 
proportional to the decay constant /q, which is known 
from heavy-quark effective theory to scale in terms of the 
heavy-quark mass mq as 1 jyjmq. In turn, the relevant 
thresholds scale as ^Sth,i ~ Tnq, so that gi ^ (sth.i)”^^"*^- 
The means by which all given coupling constants gi,meson 
are given is therefore 


5 ? = 91k 


^th,KK 


(29) 


with g^^ = 0.875 GeV^, an example given in Ref. [26] . 
Clearly, one may dispute treating the s quark as heavy, 
the accuracy of using y/Sth,i in place of mq , or indeed the 
scaling of diquark and meson cusp functions to have the 
same peak height. However, all of these issues may be 
adjusted to match one’s prejudices by including appropri¬ 
ately chosen 0(1) correction coefficients. Our prescrip¬ 
tions are designed to make direct comparisons between 
the forms as clear as possible. 

In Figs. and [^ we compare the pole-dragging ef¬ 
fectiveness parameter y obtained for diquark thresholds, 
alongside the mesonic form with the KK, DD*, and BB* 
thresholds, respectively, each with /3i = 1.0 GeV. One 
sees that the KK mesonic curve is much wider than the 
diquark curve, while the DD* and BB* mesonic curves 
are much narrower. We infer that the presence of a 6-S 
threshold is more effective in dragging a resonant pole 
for heavy-quark systems, and the two-meson threshold is 
more effective in dragging a resonant pole for light-quark 
systems. 

But several other interesting results follow from Figs.j^ 
[5] and[6j First, the maximum height of the pole-dragging 
function y is much larger in the KK case than in the 
DD* or BB* cases, for either diquark or mesonic forms. 
Even taking into account the larger value of ■s/siKi for the 



FIG. 4: Comparison of the effectiveness of pole dragging 
by cusps [y of Eq. (281] as a function of a; = Mq/ y/saD, 
from the diquark cusp (solid, green) and from the mesonic 
cusp (dashed, red) chosen to have /3* = bO GeV. Here, 
^Sth.i = 0.991 GeV, and gi,diquark = 1.348 GeV, while 
gi,meson = 0.935 GeV to give the diquark and meson cusp 
functions the same height. 



FIG. 5: As in Fig. [^ except with .^Sth.i = 3.872 GeV, and 
gi,diquark = 0.370 G^, while gi,meson = 0.474 GeV to give the 
diquark and meson cusp functions the same height. This case 
can be compared directly to that with (5 = 1 GeV in Fig. [^ 
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FIG. 6: As in Fig. H except with ^Stb,i = 10.604 GeV, and 
ffi,diquark = 0.136 GW, while Si, meson = 0.286 GeV to give the 
diquark and meson cusp functions the same height. 


in both relative and absolute terms for lighter systems 
{KK vs. DD* or BB*). We also found that the effect 
due to mesonic cusps is larger than that due to diquark 
cusps for the KK threshold, while the diquark cusp effect 
is stronger for DD* or BB* thresholds. 

This calculation is of course an idealization, in that 
only one threshold is present in each example; in real¬ 
ity, several thresholds, each with its own strength and 
contributing with different signs (due to the different 
Riemann sheets) are simultaneously present and all con¬ 
tribute to the total effect. Furthermore, this calculation 
assumes, for maximum clarity of comparison, that the 
coupling constants for the diquark and mesonic cases give 
cusp functions of the same height, and that the coupling 
constants scale with the thresholds according to expec¬ 
tations from heavy-quark effective theory. Any of these 
approximations can be relaxed in a more detailed analy¬ 
sis. 

The central conclusion, however, is that if S-S states 
exist in the spectrum of QCD, then the opening of their 
production thresholds produce measurable shifts in the 
masses of resonances, which must be taken into account 
in precisely the same way as shifts appearing due to the 
opening of meson production thresholds. The cusp effect 
appears to promise a rich source of new physical effects. 


latter cases, the absolute size of the maximal pole drag¬ 
ging is larger for the lighter cases; to wit, the numbers 
are about 200 MeV, 4.3 MeV, and 0.2 MeV, respectively. 
Some, but not all, of this decreased effectiveness can be 
attributed to the decrease of gi via heavy-quark scaling 
given by Eq. (291. Indeed, for a fixed value of Sth.i, the 
relative size of the pole-dragging effect y is found empiri¬ 
cally to scale approximately as gf ; clearly, the magnitude 
of i/sthT matters as well. In a related effect, the location 
of the peak in y travels a greater distance from thresh¬ 
old for lighter systems, and it is somewhat larger for the 
mesonic than the diquark effect, even when the cusp func¬ 
tions themselves have the same peak height. However, 
one should not conclude from these facts that the cusp 
effect is intrinsically less effective for heavy-quark sys¬ 
tems; indeed, a number of the mass splittings between 
heavy hadrons scale as l/mg, so a full analysis would 
require one to take into account values of Mq that lie 
naturally closer to y^Sth^. 


VII. DISCUSSION AND CONCLUSIONS 
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Appendix: Expressions for Unequal Masses 


In the case ^ rn 2 ,i, let us define 

m = + 7712 ,i), 

5 = - 77l2.i) , 

(52 


(A.l) 


and henceforth suppress the index i. Note that sth,i = 
Am?. The kinematical variables are, in analogue to 
Eqs. ([^ and ([^, 


We have performed the first analysis of the cusp ef¬ 
fect due to the opening of diquark-antidiquark thresh¬ 
olds by using constituent counting rules to model their 
production form factor. We directly compared our re¬ 
sults to those obtained from employing a frequently-used 
phenomenologically-based meson form factor, and found 
that the magnitude of the pole-dragging effect is greater 



and its inverse reads 

2777 ^ 

s=--^(l + e)(l + /7), (A.3) 

I — 
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written in terms of the auxiliary variable 


and that of Eq. (13) is 


h = i l- 


4e(l — p^) 


( 1 - 


(A.4) 7r(s) = 


which equals 1 in the equal-mass case e = 0. Analogous 
primed forms hold for s —> s'. 

Then the generalization of Eq. 0 becomes 


7r(s) = — P /" dp'^ 

'o 


F^s') 


TT 


h' 


( Ihp! _ ’ 

y y+Zi_ \-^h' j 


(A.5) 


dp' 


F\s') + 


4to^p'^( 1 -I- e)(l -I- h') dF'^{s') 


(l-p'2)2/i' 


X In 


2 \ /I /2n 1 + ^ 


(l-p^)-(l-p'^) 


1 + h' 


ds' 
(A.6) 
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